Gravitational self-force on a particle orbiting a Kerr black hole 
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We present a practical method for calculating the gravitational self-force, as well as the electromag- 
netic and scalar self forces, for a particle in a generic orbit around a Kerr black hole. In particular, 
we provide the values of all the regularization parameters needed for implementing the (previously 
introduced) mode-sum regularization method. We also address the gauge-regularization problem, 
as well as a few other issues involved in the calculation of gravitational radiation-reaction in Kerr 
spacetime. 



o 
o 

(N 



> 
m 
O 

(N ■ 

(N : 

o ■ 

cr 

i ' 



X 



04.70. Bw, 04.25.Nx 

The gravitational two-body problem is a long-standing 
open problem in General Relativity. In Newtonian the- 
ory the evolution of a binary system is easily analyzed, 
yielding periodic orbits with conserved energy and angu- 
lar momentum. In General Relativity, in contrast, the or- 
bits in a binary system are never periodic, and the orbital 
energy and angular momentum decrease due to emission 
of gravitational radiation. Typically, the two masses will 
gradually inspiral towards each other, until they merge. 
An interesting branch of the problem refers to the case 
where one of the components is a compact object much 
lighter than the other. In this situation, the lighter ob- 
ject may be treated as an infinitesimally small "particle" 
moving on the fixed gravitational-field background of the 
heavier object. The particle then experiences a self force 
(SF) due to interaction with its own gravitational field 
(which is distorted by the field of the massive object). 
The issue, then, becomes that of determining the gravi- 
tational SF acting on the particle. The electromagnetic 
(EM), flat-space analogue of this problem was explored 
many decades ago by Lorentz [1] and Dirac [2], who an- 
alyzed the force acting on an accelerating point charge 
due to interaction with its own EM field. To analyze the 
evolution of a strong-gravity, large mass-ratio binary sys- 
tem, the classic works by Lorentz and Dirac need to be 
generalized in two important respects: From the EM to 
the gravitational SF, and from flat to curved spacetime. 

DeWitt and Brehme [3] were first to extend the SF 
concept to curved spacetime, in the EM case. More re- 
cently, a theoretical framework has been established [4] 
for the gravitational SF in curved spacetime (cf. [5]). An 
essential feature of both the EM and gravitational SFs in 
curved spacetime (unlike the analogous flat-space case) is 
that they depend, at any given point along the particle's 
worldline, upon the entire past history of the particle. 
This may be attributed to the scattering of the (EM or 
gravitational) waves off the spacetime curvature. So far, 
this non-locality rendered actual calculations of the SF 
in curved spacetime most challenging. 

The quest for a practical SF calculation method has 



considerably intensified in recent years, with the deploy- 
ment of a new generation of gravitational-wave detec- 
tors. One of the main sources for the planned space- 
based gravitational-wave detector LISA [6] would be the 
inspiral of compact objects [white dwarfs, neutron stars, 
and stellar black holes (BHs)] into super-massive BHs, 
with a typical mass ratio of 10 3 -10 7 . Such super-massive 
BHs appear to reside in the centers of many galaxies, 
including our own [7]. In a typical scenario, the com- 
pact object spends the last year of inspiral deep inside 
the highly rclativistic region near the BH's event horizon, 
emitting some 10 5 cycles of gravitational waves [8], which 
carry a detailed information of the massive BH's strong- 
field geometry [9] . Knowledge of the self force — which, in 
turn, determines the orbital evolution — is necessary for 
predicting the emitted waveform. Such a prediction is 
also important for the design of templates which should 
significantly improve the chance of detection in the situ- 
ation of a small signal-to-noise ratio, typical to the first 
generation of detectors. 

In several, relatively simple situations the orbit's evolu- 
tion may be inferred using standard energy-momentum 
balance methods [10]. Namely, one calculates the flux 
of energy and angular momentum carried away by the 
gravitational waves, and equates this to (minus) the rate 
of change in the corresponding entities associated with 
the orbiting particle itself. This simple method yields 
the orbital evolution for all orbits around a non-spinning 
(Schwarzschild type) BH, and also for strictly equato- 
rial or circular orbits around a rotating (Kerr type) BH. 
However, in the astrophysically relevant situation one 
needs to analyze the evolution of generic orbits around 
a Kerr BH. In this generic case, energy-momentum bal- 
ance methods are insufficient, because the energy and 
(azimuthal) angular-momentum no longer uniquely de- 
termine the orbit. One therefore must calculate the local 
SF acting on the particle. The goal of this paper is to 
provide a practical method for calculating this SF, for 
generic orbits in Kerr spacetime. 

Based on the basic formulation of Refs. [4], we previ- 
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ously devised the mode-sum method, as a practical tech- 
nique for SF calculations [11,12]. In this method, one 
first calculates the (finite) contribution to the SF due 
to individual multipole modes of the particle's field (this 
is done by integrating numerically the mode-decomposed 
linearized field equations) . Then, a certain regularization 
procedure is applied to the mode sum [see Eq. (4) below] , 
which involves certain "regularization parameters" (RP). 
The latter depend upon the background spacetime and 
the specific orbit under consideration. The RP values 
were previously derived for equatorial orbits around a 
Schwarzschild BH [13-16]. Here we extend the analysis 
and give the RP values for an arbitrary geodesic orbit 
around a Kerr BH. 

Throughout this Letter we use units where c = G = 1 , 
metric signature — | — h+, and Boyer-Lindquist coordi- 
nates t, r, 9, ip. For completeness, we shall address here 
not only the gravitational SF, but also the analogous EM 
and scalar phenomena. We shall thus consider the follow- 
ing setup: A particle of mass /i, electric charge e, and/or 
scalar charge q moves along an arbitrary free-fall orbit 
(i.e., no external forces) z m (t) around a Kerr BH with 
mass M » /j,, \q\, \e\ and arbitrary spin aM. We shall 
consider the SF acting on the particle (to the leading or- 
der oc , e 2 , q 2 ) at an arbitrary point Zq = (io, ^o, 9q, ¥>o) 
along its trajectory. The particle's equation of motion, 
including the SF effect, reads 

/xu^u" = Ff f , (1) 

where u 13 = dz 13 /dr is the four- velocity, a semicolon de- 
notes covariant differentiation with respect to the back- 
ground (Kerr) geometry, and F| elf is the SF. 

Let x denote a spacetime point in the neighbor- 
hood of Zq. We denote the particle's actual perturba- 
tion field (the "full field") by h a p{x), A a (x), or $(a;), 
representing, correspondingly, the metric perturbation 
(oc ^i), EM vector potential (oc e), and scalar field 
(oc q). These fields are, in principle, solutions to the lin- 
earized Einstein/Maxwell/Klein-Gordon equations, with 
the particle serving as a source term. To these fields there 
correspond "full-force" fields, given by [14,15] 

{q$,p, (scalar), 
eKfiAp„, (EM), (2) 
HK^hp^s, (grav.), 

where the tensors k are 

K a ^(x)=g a ^u> 3 -g aP u<, 

K afh *(x) = -g a(3 u<u 5 + ^g^iiPu 1 - ^u^u^u 5 . (3) 

Here g a @ is the background metric at x and u a is the 
four-velocity parallelly propagated from the worldline to 
x along a short geodesic normal to the worldline. 

In the mode-sum method, the SF is then calculated 
through [11] 



K eli = E [(tX)* A t L - B » - c »i L \ D »- ( 4 ) 

1=0 

Here L = I + 1/2, and (Fffi 11 )* is the I multipole mode 
of Fr 1 (x), summed over all azimuthal numbers m (for 
a given I) and evaluated at (to,r — ► r^, 6o, </?o)- The 
full modes are generally discontinues across r = r , 
hence the label ± [11]. The RP, i.e., the quantities A^, 
Bp, Cp, Dp, are to be derived analytically (by analyzing 
the local singular behavior of the perturbation field near 
the particle) . The first three of these parameters capture 
the asymptotic form of the full- force modes at large I, 
making the sum in Eq. (4) convergent. The parameter 
Dp is a certain residual quantity that arises in the sum- 
mation over I — for a more precise definition of the RP 
sec [11,13]. 

In deriving the RP in Kerr spacetime, we basically fol- 
lowed the method used in the Schwarzschild case [13], 
as described in much detail in Rcfs. [14,15]. Namely, 
we considered the direct part of the full- force field (i.e., 
the part associated with waves directly propagating along 
the light cone), from which all RP are obtained by apply- 
ing the multipole decomposition and appropriately tak- 
ing the limit x — > z$. The main technical challenge in 
extending our analysis to the Kerr case lies in the fact 
that one can no longer exploit spherical symmetry by 
choosing a coordinate system in which the orbit is con- 
fined to the equatorial plain. Rather, one must carry out 
the calculation at an arbitrary value of 6q. This turns 
out to render the entire analysis significantly more com- 
plicated. For lack of space, we skip the derivation of 
these RP values — the detailed derivation will be given 
elsewhere [17]. In what follows we merely present the re- 
sults, i.e., the values of all RP for an arbitrary geodesic 
in Kerr spacetime. 

To write down the RP in a convenient unified form, we 
introduce a "generalized" charge q s , where q s — q, e or 
/i, for the scalar (s = 0), EM (s = 1), and gravitational 
(s = 2) cases, respectively. All RP are then independent 
of s, apart from an overall factor (—l) s q 2 . For A^, C M , 
and Dp we find 

Cp = Dp = 0, (5) 

and 




Af = -(u r /u t )Af, A± = A±=0, (6) 

where hereafter g a p = g a/3 (z a ) and V = 1 + u 2 e /g g0 + 
u %/9f<P- The expression for Bp takes a slightly more 
complicated form, 

Bp = (-l) s q 2 s {2^)- 1 Pp abcd I abcd . (7) 
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Hereafter, roman indices (a, b, c, ...) run over the two 
Boyer-Lindquist angular coordinates 9, ip. The coeffi- 
cients P^abcd are given by 

Pfiabcd = {<^P[iaPbe P[iePab)C c d 

+ \ [SP^dPabc - (2P^b + Pab^)Pcd] , (8) 

where 

Paf) = 9a/3 +U a U(3, 

P a f3j = (u\u 7 r*p + g a p, 7 /2) (9) 
(r^ denoting the connection coefficients at z ), and 

C^ = isin0ocos0 o , = C* e = -\ cot0 o , (10) 

with all other coefficients C e cd vanishing. Finally, the 
quantities I abcd are 

I abcd = (sm6 )- N / 2 "G( 7 )- 5/2 (sin 7 ) Ar (cos7) 4 - A, d7, 
Jo 

where TV = N(abcd) is the number of times the index if 
occurs in the combination abed (namely, N — 5® + 6^ + 
5° v + S d ) and 

G(7) = Pq,q, sin 2 7 + 2Pg^ sin7COS7 + Pee cos 2 7, (11) 
where 

P^0 = P vv / sin 2 60, P&0 = Pe v / sin (12) 

The integrals I abcd can be expressed in terms of complete 
Elliptic integrals. These relations are given explicitly in 
Ref. [18]. 

In the special case of an equatorial orbit in 
Schwarzschild spacetime, the above RP values reduce to 
those obtained previously using various methods [11,13]. 

In the rest of this paper we discuss several important is- 
sues that arise when implementing the mode-sum method 
in Kerr spacetime. 

First, to avoid confusion, we emphasize that the quan- 
tities Ff^ n appearing in Eq. (4) refer (in all three cases 
s = 0,1,2) to the decomposition of the full-force field 
components F^ u11 in the standard scalar spherical har- 
monics. Practically, in the Kerr case one separates 
the field equations in the frequency domain, using (spin 
weighted) spheroidal harmonics. Then the full-force field, 
too, is obtained in terms of the spheroidal harmonics. 
Thus, in order to implement the mode-sum scheme (at 
least in its present formulation), one needs to decom- 
pose the contributions to the full force from the vari- 
ous spheroidal-harmonic modes, into scalar spherical har- 
monics: For a given I (and given azimuthal number m 
and temporal frequency to) one collects the contributions 
to the spherical- harmonic mode I, m coming from the 



full-force fields associated with the spheroidal-harmonic 
modes l',m,u. Then one sums over I' (as well as m, w) 
to obtain Fjfi 11 . The decomposition of the spheroidal har- 
monics in spherical harmonics is described in Ref. [19]. 

Our second remark concerns the construction of the 
full-force field in the EM and gravitational cases. This 
field involves an extension of u a off the worldline, cf. Eq. 
(3). Here we adopted the most natural extension u a , 
namely, the one obtained by parallelly propagating u a 
from the worldline to x. As discussed in [15], it is pos- 
sible to use any other (sufficiently regular) extension of 
u a . Note, however, that changing the extension (which 
obviously affects the quantities F^ u ) will generally affect 
the value of the parameter (though it turns out that 
the other RP are unaffected). Nevertheless, if two exten- 
sions share the same values of angular derivatives u a b at 
the worldline, they admit the same [15]. 

The parallelly-propagated extension is an elegant one, 
and, furthermore, it allows the simplest derivation of B^. 
However, this extension is rather inconvenient for calcu- 
lating the full-force quantities (the Legendre decom- 
position requires the extension of u a far away from the 
worldline, on the entire two-sphere r, t=const; It is hard 
to calculate u a explicitly off the worldline, and it is even 
harder to Legendre-decompose it. Furthermore, there is 
no guarantee that u a would be globally well defined) . To 
overcome this difficulty, we propose the following strat- 
egy: Compose any global extension u a off the worldline, 
with the only demand that u a b — u a b (and, obviously, 
u a = u a ) at the worldline itself. It is easy to construct 
such simple global extensions. Then simply use this u a 
instead of u a in Eq. (3). The above values of all RP are 
unaffected by this change of extension. 

There is another technical issue that arises in the EM 
and gravitational cases: The only known scheme for sep- 
arating the EM and (linearized) gravitational field equa- 
tions in Kerr spacetime is the Teukolsky formalism [20]. 
By solving the Teukolsky equation, one obtains certain 
components of the Maxwell or Weyl tensors, respectively. 
However, as is obvious from Eq. (2), the calculation of 
F^ n (and hence of F^f 1 and the SF) requires the knowl- 
edge of the basic perturbation fields, h a p(x) or A a {x). 
A method for reconstructing h a p and A a for a point 
particle, out of the Teukolsky variables, based on the 
Chrzanowski-Wald formalism [21,22], was recently pro- 
vided by Ori [23] (see also [24] for the Schwarzschild case) . 
More recently, this construction was extended to the non- 
radiative modes I = 0, 1, as well as to all higher-Z modes 
with u = [25]. 

Finally, we address the crucial issue of the gauge- 
dependence of the gravitational SF, and the related prob- 
lem of gauge regularization in Kerr spacetime. The gravi- 
tational SF (unlike the scalar and EM forces) turns out to 
be gauge-dependent [26] . The basic formalism [4] is given 
in terms of the harmonic gauge. The transformation of 
the SF to any other gauge is given in Ref. [26] . This trans- 
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formation, however, is restricted to regular gauges, i.e., 
gauges related to the harmonic gauge via a displacement 
vector £^ that is continuous on the particle's worldline. 
(If the gauge is irregular, then the SF generally becomes 
ill-defined.) For any regular gauge, the SF may, in prin- 
ciple, be calculated through Eq. (4) — provided that one 
uses the quantities Ffjf 1 associated with this same gauge. 
The RP are gauge- invariant in this context [26] . 

Although the momentary SF is gauge-dependent, the 
long term radiation- reaction evolution of the orbit, as 
expressed, e.g., by the drift of the constants of motion, 
is gauge- invariant. The radiation-reaction evolution can 
therefore be calculated in any regular gauge. 

Chrzanowski's construction yields the metric perturba- 
tions in the so-called radiation gauge [21]. Unfortunately, 
this gauge is not regular, hence the SF does not have a 
well-defined value in the radiation gauge. To overcome 
this difficulty, we previously proposed [26] to calculate 
the SF in an intermediate gauge, namely, a regular gauge 
which (unlike the harmonic gauge) can be obtained from 
the radiation gauge via a simple, explicit gauge trans- 
formation. A convenient such intermediate gauge was 
recently constructed explicitly in Ref. [25]. 

When implementing the mode-sum scheme (4) to the 
SF in the intermediate gauge, one should, in princi- 
ple, substitute the quantities _p^ u ( INT ) associated with 
this same gauge. Alternatively, one may still use the 
radiation-gauge quantities p^ n ( RAD ) m gq but in 
this case one must subtract an extra term, which we de- 
note 5 Dp, to compensate for the differences SF^ U in the 
values of F^f 1 in the two gauges. (The derivation of SF^ 11 
from £^ is outlined at the end of Ref. [26]; SD^ is the sum 
of SFjj 11 over I.) Thus, when applied to the SF in the in- 
termediate gauge, Eq. (4) may be recast as 



^sclf(INT) _ 

1=0 



fa 



full(RAD) 

fj.1 



B„ 



5D„ 



(recall C M = = 0). The RP A^,B^ are those given 
in Eqs. (6) and (7) (as mentioned above, the RP values 
are gauge- invariant). The term SD^ was recently cal- 
culated for a generic orbit in Kerr, and its expression 
(whose structure and complexity resemble those of 
above) will be given elsewhere [25]. This provides the 
solution to the gauge-regularization problem. The quan- 
tity ^eif(iNT) can now be usec i to analyze the (gauge- 
independent) long-term radiation-reaction evolution, for 
any orbit in Kerr spacetime. 
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